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New constraints on dark energy from the observed growth 
of the most X-ray luminous galaxy clusters 
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ABSTRACT 

We present constraints on the mean matter density, Q m , the normalization of the 
density fluctuation power spectrum, as, and the dark-energy equation-of-state param- 
eter, w, obtained from measurements of the X-ray luminosity function of the largest 
known galaxy clusters at redshifts z < 0.7, as compiled in the Massive Cluster Survey 
(MACS) and the local BCS and REFLEX galaxy cluster samples. Our analysis employs 
an observed mass-luminosity relation, calibrated by hydrodynamical simulations, in- 
cluding corrections for non-thermal pressure support and accounting for the presence 
of intrinsic scatter. Conservative allowances for all known systematic uncertainties are 
included, as are standard priors on the Hubble constant and mean baryon density. 
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model, and fi m = 0.24Z 07 , ^ 8 = 0.851^ and w = -1.4 

model (marginalized 68 per cent confidence intervals) . Our findings constitute the first 
determination of the dark-energy equation of state from measurements of the growth 
of cosmic structure in galaxy clusters, and the consistency of our result with w = — 1 
lends additional support to the cosmological-constant model. Future work improving 
our understanding of rcdshift evolution and observational biases affecting the mass- 
X-ray luminosity relation have the potential to significantly tighten these constraints. 
Our results are consistent with those from recent analyses of type la supernovae, 
cosmic microwave background anisotropies, the X-ray gas mass fraction of relaxed 
galaxy clusters, baryon acoustic oscillations and cosmic shear. Combining the new X- 
ray luminosity function data with current supernova, cosmic microwave background 
and cluster gas fraction data yields the improved constraints Q m = 0.269 ± 0.016, 
cr 8 = 0.82 ± 0.03 and w = -1.02 ± 0.06. 
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1 INTRODUCTION 

In the hierarchical collapse scenario for structure forma- 
tion in the universe, the number density of collapsed ob- 
jects as a function of mass and cosmic time is a sensitive 
probe of cosmology. The galaxy clusters that occupy the 
high-mass tail of this population provide a powerful and 
relatively clean tool for cosmology, since their growth is 
predominantly determined by linear gravitational processes. 
In the past, the local population of galaxy clusters has 
been used to jointly constrain the average matter density 
of the universe a nd the amplitude of perturbations in the 
density field (e.g. iReiprich fc Bohringerl [20021: ISeli 
IViana et all 120021 : lAUen et all 120031 : IPierpaoli et 



Schuecker et~aH 120031; IVoevodkin fc Vikhlininl liool ; iDahk 
20061 ; I Rozo et al. I l2007l ~ Pushing observations to higher 



redshift br eaks the degeneracy be t ween those two param- 
eters (e.g. iDonahue fc VoilJ Il999l; lEke etafl ll99Sl ; Irlenrvl 
l200d ; iBorgani et all l200ll ; IVikhlinin et al.l l2003h . and al- 
lows propertie s of d ark energy to be p robed as wel l (e.g . 



Haiman et all l200ll; iLevine et alj l2002t IWeller et all [20021 ; 



Maiumdar fc Mohrll2003l , 12004 ; lHenry| |2004) . 
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Investigations of this type require sky surveys with well 
understood selection functions to find clusters, as well as a 
relation linking cluster mass with an observable. A success- 
ful solution to the former requirement has been to identify 
clusters by the X-ray emission produced by hot intraclus- 
ter gas, not ably using data from the ROSAT All-Sky Sur- 
vey (RASS; lTrumpeifll993h. The ROSAT B rightest Clus- 
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Ebeling et al"1 [l998, 2000() and ROSAT- 
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ESQ Flux Limited X-ray sample (REFLEXi lBohringer et al.1 
120041 ') together cover approximately two-thirds of the sky out 
to redshift z ~ 0.3 and contain more than 750 clusters. The 
Massive Cluster Survey (MACS; lEbeling et al.ll200ll . I2007T I 
- which at this writing contains 126 clusters and covers 55 
per cent of the sky - extends these data to z ~ 0.7. 

The most straightforward mass-observable relation to 
complement these X-ray flux-limited surveys is the mass- 
X-ray luminosity relation. For sufficiently massive (hot) ob- 
jects at the relevant redshifts, the conversion from X-ray 
flux to luminosity is approximately independent of tem- 
perature, in which case the luminosities can be estimated 
directly from the survey flux and the selection function is 
identical to the requirement of detection. In a flux-complete 
survey further restricted to high luminosities, every cluster 
should thus be usable in the analysis, without the need for 
additional observations other than those required to cali- 
brate the mass-luminosity relation. A disadvantage is that 
there is a large scatter in cluster luminosities at fixed mass; 
however, sufficient data allow this scatter to be quantified 
empirically. Alternative approaches use cluster temperatur e 
l|Henrvl200d ; ISeliakl2002];|Pierpaoli et al.ll2003l ; lHenry|2004h , 
gas fraction JVoevodkin fc Vikhlininll2004l ') or Yx parameter 
l|Kravtsov et aij2006h to achieve tighter mass-observable re- 
lations at the expense of reducing the size of the samples 
available for analysis. The need to quantify the selection 
function in terms of both X-ray flux and a second observ- 
able additionally complicates these efforts. 

In this paper, we use the observed X-ray luminosity 
function to investigate two cosmological scenarios, assuming 
a spatially flat metric in both cases: the first includes dark 
energy in the form of a cosmological constant (ACDM) ; the 
second has dark energy with a constant equation-of-state 
parameter, w (wCDM). In the latter case, we account for 
the evolution of density perturbations in the dark-energy 
fluid, assuming that the dark-energy sound speed is equal 
to the speed of light. For each model, our results are in 
good agreement with findings from independent cosmologi- 
cal data sets, notably type la supernovae (SNIa), the cosmic 
microwave background (CMB), the X-ray gas mass fraction 
of galaxy clusters (/gas); an d measurements of cosmic shear. 

The theoretical background for this work is reviewed in 
Section [21 Section [3] details the data used to constrain the 
mass-luminosity relation and the cluster samples used to 
measure the X-ray luminosity function. The analysis proce- 
dure is described in Section [J] and the cosmological results 
are presented in Section [5] The importance of various sys- 
tematic effects is discussed in Section [S] 

Unless otherwise noted, masses and luminosities quoted 
in this paper or shown in figures are computed with respect 
to a spatially flat ACDM reference cosmology with Hubble 
constant h = H /W0 kms" 1 Mpc' 1 = 0.7 and fi m = 0.3. 
Luminosities and fluxes refer specifically to the 0.1-2.4 keV 
energy band in the source and observer rest frames, respec- 
tively. We will consistently use the notation L to denote to 
the true luminosity of a cluster and L to denote the luminos- 
ity inferred from observation. We will also write, for exam- 
ple, Qm to refer to the present day matter density in units 
of the critical density, whereas flm(z) is the same quantity 
at redshift z. 



2 THEORY 

The variance of the linearly evolved density field, smoothed 
by a spherical top-hat window of comoving radius R, enclos- 
ing mass M — 47rp m J? 3 /3, is 



a(M,z) 



D\z) 
2tt 2 



k 2 P{k)\W M {k)\ 2 dk. 



(1) 



Here p m is the mean comoving matter density of the uni- 
verse, P(k) is the linear power spectrum at redshift zero, 
Whiik) is the Fourier transform of the window function, and 
D(z) = a${z) I a&(0) is the growth factor of linear perturba- 
tions at scales of 8h~ x Mpc, normalized to unity at z = 0. We 
evaluate the transfer function, T(k), and pow er spectrum, 
P(k) oc k ns T 2 (k), using the CAMB package of iLewis et~ail 
(I200C " 



Ijenkins et al.l (|200ll . hereafter J01) and lEvrard et al.l 
l|2002h have shown that the predicted mass function of 
galaxy clusters of mass M at redshift z can be written in 
terms of a "universal" function of a~ 1 (M, z), 



f{? ) = ~ 



M dn(M, z) 



p m diner 1 
which can be fit by a simple form, 
/(a" 1 ) =Aexp(-| In a' 1 + B\<) , 



(2) 



(3) 



for cosmological-constant models. It has since been ver- 
ified that this fitting function is approximately univer- 
sal (within ~ 20 per cent) among m odels with constant 
w ^ — 1 and some evol v ing- w model s dKlvpin et al.l 120031; 
Linder fc Jenkins! 120031 ; lLokas et all 120041 : iKuhlen et al l 
2005h . We adopt the values A = 0.316, B = 0.67 and 



e = 3.82 from J01, determined using a spherical overden- 
sity group finder at an overdensity of 324 with respect to 
the mean matter density (324m). The number density per 
unit mass of galaxy clusters of mass M at redshift z is then 



dn(M, 



pm din a' 



(4) 



dM M dM 

Following Mor andi et al.l (|2007T ). we describe the rela- 
tionship between mass and X-ray luminosity for massive 
clusters as self-similar (e.g. iBrvan fc Normanlll998l ). mod- 
ified by an additional redshift-dependent factor 



E(z)M A = M 



E{z) 



(1 + 2)' 



(•») 



where E(z) = H(z)/Hq and Ma is the cluster mass de- 
fined at an overdensity of A with respect to the critical 
density (i.e. where the edge of the cluster is defined such 
that the mean density enclosed is a fixed factor A times 
the critical density). Our model also includes a log-normal 
scatter of width rj(z) in luminosity for a given mass. There 
is good theoretical reason to expect that this scatter may 
evolve with redshift, since some of the variability in cluster 
luminosities can be attributed to the develo pment of cool 
cores as well as the influ ence of merger events (|P'Hara et al.l 
120061 : Ichen et al1l2007l ). Since current data are insufficient 
to directly detect ev olution in the mass-luminosity scatter 
ijO'Hara et al.ll2~007l and references therein), we allow for a 
generic, linear evolution in the scatter: 

1 http://www.camb.info/ 
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r/(z) = rjo (1 + Tjiz) . 



(6) 



The form of the self-similar evolution appearing in 
Equation [5] is appropriate for constraining the mass- 
luminosity relation at fixed overdensity with respect to the 
critical density. We have chosen to constrain the relation 
at an overdensity of A = 500 in order to m atch our data 
(Secti on 13. l.lfl . and because the simulations of lEvrard et alJ 
(2008) indicate that the virial relations have lower scatter 
at this overdensity compared with higher overdensities. The 
fact that the mass function and mass-luminosity relation 
use different definitions of cluster mass requires us to convert 
masses from one overdensity to another during the analysis 
[jWhitefeOoj j; details of car rying out this conversion, assum - 
ing a spherically svmmetric lNavarro. Frenk. fc White (|l997T ) 
density profile, are reviewed in iHu fc Kravtsovl(|2003l ). 

Equations [4}{6] combined with a suitable stochastic 
model linking the true luminosity of a cluster with its ob- 
served flux, provide the means to predict the galaxy clus- 
ter luminosity function based on a set of model parameters 
(Section 0}. 



3 OBSERVATIONS 

3.1 Mass— luminosity relation 

3.1.1 X-ray data 

We have determined the mass-luminosity relation of 
galaxy clusters at lo w re dshift using the sample of 
iReiprich fc Bohringerl (|2002l , hereafter RB02), shown in 
Figure [T] In these data, the luminosities were measured 
from a combination of pointed ROSAT PSPC and RASS 
data. The masses within rsoo (M500) were determined 
by fitting the surface brightness profile with a /3-model 
|Cavaliere fc Fusco-Femianolll978T ) and applying the hydro- 
static equation, assuming that the intracluster gas is isother- 
mal. Here rsoo is the radius within which the mean density 
is 500 times the critical density, p c (not the universal mean 
matter density, p m ). 

The assumption of hydrostatic equilibrium intro- 
duces a bias into the derived masses due to the pres- 
ence of non-thermal support and asphericity in the ob- 
served clusters. Thi s bias h as been observed i n simulations 
jFaltenbacher et all 120051 ; iRasia et all 120061 ; iNagai et ail 
12007ft and in comparisons of masses derived from X-ray and 
weak gravitational len sing observations of the same targets 
l|Mahdavi et al.ll20Qg| ). alt hough the number of objects in 
these studies remains low.[R asia et al.1 l|2006l) find a bias of 
—34 ± 14 p e r cen t at rsoo from simulations of 5 clusters. 
INagai et al.1 (120071 ) find a smaller bias, —25 ± 16 per cent, 
from 16 simulated clusters at z = 0. Note that there are 
differences in the methodology of the two groups; Rasia et 
al. use an isothermal /3-model to determine the X-ray mass 
from mock images, while Nagai et al. use a more sophis- 
ticated X-ray analysis and additionally con sider the obser- 
vation al bias in the determination of rsoo- M ahdavi et al.l 
(2008) have compared masses determined from X-rays and 
lensing for 18 clusters and find a bias of —22 ± 9 per cent 
at rsoo, assuming the lensing masses are accurate. We adopt 
the Nagai et al. values, —25 ± 16 per cent, as the nomi- 
nal values of the mean bias, b, and scatter in bias between 




Figure 1. Mass— luminosity data of RB02 and the best- fitting 
relation (Equation . The quantities Y and X\ are defined in 
Equation [8] and relate to log-mass and log-luminosity, respec- 
tively. As noted in the text, objects with X\ < —0.5 or z > 0.11 
are excluded from the analysis. All masses and luminosities are 
computed for our reference cosmology. 



clusters, Sb- In order to account for the spread among the 
various results for b and the small number of simulated clus- 
ters used to estimate s^, we assign a 20 per cent systematic 
uncertainty to both of these quantities (see below). 
We fit the log-linear model (see Equation [5} 



Y = q + /3Xi+ 7 X 2 , 



(7) 



where 

y = log 

Xi = log 



E(z)M 500 



M 



© 



L 



• !:) l ^(z)10 44 ergs- 1 ' ' 
•Vi> - • ,0+2): 

, Mo 

" = < <!! -'- 



(8) 



The full RB02 data set contains some small groups which 
can be influential in the estimation of the slope, f3. Since we 
only require the mass-luminosity relation to be calibrated 
approximately one decade below the luminosity cut used for 
the X-ray luminosity function data (Section l3.2p . we fit only 
the data with Xi > —0.5 in our reference cosmology. This 
choice mitigates the possibility of obtaining biased results if 
slope of the mass-luminosity relation is different for massive 
clusters compared with smaller groups. 

The process of fitting the model in Equation [7] is com- 
plicated by the presence of Malmquist bias. Close to the flux 
limit for selection, any X-ray selected sample will preferen- 
tially include the most luminous sources for a given mass. 
This results in a steepening of the derived mass-luminosity 
relation and a bias in the inferred intrinsic scatter in lu- 
minosity for a given mass. (The inferred intrinsic disper- 
sion may be artificially reduced or increased, depending on 
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Figure 2. Luminosity-redshift distribution of the RB02 data. 
The quantity X\ is defined in Equation \E\ and is computed in our 
reference cosmology. The data from RB02 which are excluded 
from the analysis (see text), are shown as blue triangles. The 
remaining clusters have an extent of roughly one decade in lumi- 
nosity at all redshifts. 



the distribution of data with respect to the flux limit.) The 
use of the extended sample of RB02, rather than only their 
flux-limited HIFLUGCS sample, partially mitigates this ef- 
fect by softening the flux limit. Furthermore, we eliminate 
the 6 clusters at redshifts > 0.11 which all lie close to the 
flux limit; as Figure [2] shows, the extent in luminosity of the 
remaining 86 clusters is roughly a decade at all redshifts. 
As our estimate (see below) of the intrinsic scatter in Xi 
is roughly 0.12 (i.e. 0.12 decades in luminosity), and the 
Poisson measurement error is much smaller, we expect the 
effects of Malmquist bias on our estimation of the power-law 
slope and intrinsic dispersion to be minimal. 

A second consequence of Malmquist bias is a strong 
apparent, but not necessarily physical, correlation between 
luminosity and redshift due to the fact that the flux limit 
corresponds to higher luminosities at higher redshifts. (This 
is evident in Figure [2]) Within redshift 0.11, we expect this 
false signal to be much greater than any real evolution; we 
therefore fix 7 = when fitting the model. 

Although methods exist for performing linear regression 
on data wi th bivariate heteroscedast ic errors and intrinsic 
scatter (e.g. lAkritas fc Bershadvll996l ). such methods do not 
provide a simple goodness-of-fit measure that can be associ- 
ated with arbitrary values of the fit parameters. In contrast, 
the \ 2 statistic is an easily calculated measure of goodness- 
of-fit, but is biased by the presence of intrinsic scatter. We 
compromise by using a modified \ 2 statistic which accounts 
for intrinsic scatter by introducing an additional dispersion 
term, 



-2 
X 



E 



(9) 



where ey,j is the measurement error on Yj and 8 is the ad- 
ditional dispersion. This statistic is defined in terms of Y 




Figure 3. 68.3 and 95.4 per cent confidence regions for the mass- 
luminosity parameters a and (3 from the RB02 data (solid lines). 
Also shown is the best fit and 68.3 per cent confidence inter- 
val (square with error bars) on a with fixed j3 obtained from 
the Dahlc (2006) data (Section [5X2) . The cosmology is fixed at 
our reference. The dashed lines are the confidence regions that 
would have been obtained from the RB02 data if the masses were 
not corrected for bias due to the assumption of hydrostatic equi- 
librium. The difference in the precision of the a constraints be- 
tween the corrected and uncorrected RB02 results comes from our 
marginalization over the uncertainty in the size of the correction. 



given Xi because the luminosity measurement errors are 
negligible in comparison to the mass measurement errors 
(Figure [l]). The value of 8 is chosen iteratively, adjusting it 
such that the best fit, found by minimizing \ 2 using least- 
squares methods, has % 2 equal to the median of the chi- 
square distribution with v = — 2 (86 data points minus 
2 free parameters) degrees of freedom (i.e. \ 2 ' /v w 1). The 
X 2 statistic thus does not measure absolute goodness-of-fit, 
but goodness relative to the best fit. Note that in the cos- 
mological analysis we use the fixed value 8 = 0.014 from our 
reference cosmology (see Section [4.11 for details). 

A common log-normal intrinsic scatter in luminosity for 
a given mass can be estimated from the scatter in the data 
about the best fit, with measurement errors subtracted in 
quadrature. Since these data are predominantly at very low 
redshift, we use this procedure to estimate the z = value 
of the intrinsic dispersion (see Equation [6j , 



- 2 

VO = 



Y, 



OLf 



0s 



Al,-, 



2 

£ X 1 ,j 



01 



,(10) 



where Uj and (3f are the parameters describing the best 
fit. Under assumptions of normal measurement errors and 
intrinsic scatter, the quantity rjefv/rj 2 . is drawn from a chi- 
square distribution with v degrees of freedom. This likeli- 
hood for the measured value rjo can then be used to constrain 
the model parameter rjo- 

The scatter in mass bias between clusters, St,, is prop- 
agated into the mass error bars, ey, and thus impacts di- 
rectly the estimated intrinsic scatter, rjo (Equation llOp . This 
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results in a degeneracy between s& and T70 ; however, Sb is not 
coupled to the other fit parameters, and so we can take its 
systematic uncertainty into account simply by increasing the 
width of the likelihood function for 770. We therefore replace 
the chi-square likelihood discussed above by the gamma dis- 
tribution (of which the chi-square distribution is a special 
case) . Fixing the location of the distribution's peak and aug- 
menting its variance by a factor of F, the shape and scale 
parameters of the gamma distribution, k and 6, can be de- 
termined by solving the system of equations 

(k -1)0 = v-2, 

k6 2 = {2v)F. (11) 

In practice, the value F — 14 approximately corresponds to 
the desired 20 per cent uncertainty in st- Finally, we account 
for a 20 per cent systematic uncertainty in the mean bias, 
b, by straightforwardly marginalizing over an appropriate 
prior. In detail, we parametrize the mean bias through B = 
— log 10 (l — b)= log 10 (M truo /M cst ), since this quantity is 
added to the log-masses, Y. 

Confidence regions for a and /3 obtained with cosmo- 
logical parameters fixed at our reference and using uniform 
priors on a, /3 and rjo are displayed in Figure [3] (solid lines). 
The best-fitting values and marginal 68.3 per cent confi- 
dence intervals are a = 14.49 ± 0.04, f3 = 0.671^ and 
770 = 0.121q qj. A quantile-quantile plot of the residuals in 
X\ for the best fit (Figure [4j confirms that the distribution 
of luminosities for a given mass is reasonably approximated 
by the log-normal distribution 

3.1.2 Weak lensing data 

In order to verify the appropriateness of the hydrostatic bias 
correction applied to the RB02 masses, we com p are th e re- 
sults of Section |3~TT1 with the data of iDahld (120061 ). for 
which masses were measured using weak gravitational lens- 
ing. These masses, with associated luminosities taken from 
the BCS and extended BCS (eBCS) catalogs, are shown in 
Figure [5] (black points). In contrast to the RB02 data, the 
Malmquist bias due to the (e)BCS flux limit is clearly ev- 
ident; in particular, the least massive clusters in this data 
set are very likely to represent the upper tail of the distri- 
bution of luminosities for given mass. Although the simple 
methods described above are unsuitable for fitting the mass- 
luminosity relation when the data are subject to significant 
Malmquist bias, the intersection of these two data sets near 
the highest observed masses indicates qualitatively that they 
are compatible. 

W e can quantify this observation by using the IDahld 
(2006) data to fit for the normalization, a, by minimizing 
X 2 , while leaving the slope, (3, fixed at the best-fitting value 
from Section 13.1.11 The resulting one-dimensional 68.3 per 
cent confidence interval is a = 14.461q'qs (Figure in 
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Figure 4. Quantile-quantile plot comparing the residuals in X\ 
of the RB02 data for the best fitting relation to the normal dis- 
tribution. The straightness of the distribution indicates that the 
scatter is well approximated as normal in X\ (log-normal in lu- 
minosity). (See footnote 2 .) Inset: the histogram of X\ residuals 
is directly compared with the normal distribution. 




I I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 

-0.5 0.5 1 1.5 

Figure 5. The mass-luminosity data of lDahlj l|2006t ) (black) are 
compared with those of RB02 (blue; see additional comments of 
Figure [TJ. The effect of the (e)BCS flux limit on the Dahle data 
(X\ ~ 0.75) is clearly evident. At lower masses (Y), the data set 
is progressively more likely to contain only the most luminous 
clusters for that mass, resulting in an apparent steepening of the 
relation. 



2 The quantile function for a probability distribution is the in- 
verse of the cumulative distribution function (CDF). Figure [4] 
shows each residual ri against the quantity F" 1 (F^(ri)), where 
Fjv is the standard normal CDF and Fe is the empirical CDF 
of the residuals. If Fe is approximately normal, the points will 
be distributed roughly along a line. Significant curvature would 
indicate a poor match in distribution. 



agreement with the corrected RB02 data. In contrast, the 
lensing constraint is incompatible at the 68 per cent con- 
fidence level with the confidence regions that would have 
been obtained from the X-ray data without applying the 
bias correction (dashed lines in Figure [3]). 
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3.2 X-ray luminosity function 

We use three flux-li mited surveys in our a nalysis: the BCS 
llEbeling et al.lll998f ) and REFLEX sample i|Bohringer et all 
l2004h at low reds hifts (z < 0.3), and the MACS 
labeling et al J l200ll ) at 0.3 < z < 0.7. (We restricted the 
REFLEX sample to the southern hemisphere so that its 
coverage on the sky would not overlap the BCS.) As dis- 
cussed in Section 14.21 a proper accounting of the intrinsic 
scatter in the mass-luminosity relation involves convolv- 
ing over all possible masses when evaluating the likelihood 
of a set of cosmological parameters; it is therefore neces- 
sary for the mass-luminosity relation to be well calibrated 
at luminosities significantly (at least an order of magni- 
tude) below those allowed in the surveys. We thus also re- 
strict the analysis to clusters with large inferred luminosities 
L > 2.55 x 10 44 hf 2 ergs -1 . 

The completeness of the R EFL EX sample was i nvest i- 
gated bv lBohringer et al.l (|200lh and lSchuecker et all i|200ll ). 
and is thought to be > 90 per cent at a flux limit of 
3.0 x 10 -12 ergs -1 cm -2 . T he BCS completeness as a func- 
tion of flux is quantified in lEbeling et al.l (|l998l ); we use a 
flux limit of 4.4 x 10 -12 ergs -1 cm -2 where the BCS com- 
pleteness is 90 per cent. Most of the incompleteness in the 
BCS is due to the inefficient extended-source detection of 
the RASS II algorithm, which is most severe at very low red- 
shifts. One consequence of our high luminosity cut is that the 
sample contains mostly higher redshift objects (z w 0.21) 
and not the large number of low redshift, low luminosity ob- 
jects that would be included in a strictly flux-limited sam- 
ple. The reported completeness of the BCS is thus a sig- 
nificant underestimate for the sub-sample of very luminous 
clusters used in our analysis. We have repeated the analy- 
sis of the BCS data alone (see Section I5.1jl using a higher 
flux limit (5 x 10 -12 ergs -1 cm -2 ) where the reported BCS 
incompleteness is negligible, and obtain a similar result. We 
conclude that significant incompleteness is not present in 
the BCS, given our selection criteria. The similar numbers of 
clusters (78 and 80, respectively) in the BCS and REFLEX 
samples satisfying the BCS flux limit and our luminosity cut 
also support this conclusion. 

Unlike the BCS and REFLEX, for which extended- 
source fluxes were measured using the Voronoi tessella- 
tion a nd percolation (|Ebelindll993l; lEbeling fc Wiedenmaiml 
1 19931 ) and growth curve analysis (|Bohringer et al.ll2000l ) al- 
gorithms, respectively, reported MACS fluxes are measured 
within a fixed angular size aperture (5 arcmin in most 
cases). A r edshift-dependent co rrection for missing flux, de- 
scribed in E beling et al.l l|200lf) . is required when convert- 
ing this aperture flux to total flux, both when determining 
luminosities for the detected clusters themselves and when 
computing the apert ure flux-dependent survey sky coverage 
l|Ebeling et al.ll2007l ). We adopt an aperture flux limit of 
2 x 10 -12 ergs -1 cm -2 , for w hich the completeness is > 90 
per cent (|Ebeling et alj|200ll l. 

Using the flux and redshift limits and the luminosity 
cut discussed above, the BCS, REFLEX and MACS sam- 
ples respectively contribute 78, 130 and 34 clusters to our 
sample. The luminosity-redshift distribution of these data 
is displayed in Figure [IJ] 

In order to adequately take the effects of Eddington bias 
into account when predicting the number of clusters above 
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Figure 6. Luminosity-redshift distribution of clusters in the BCS 
(blue, filled circles), REFLEX (black, open circles) and MACS 
(red triangles) which are above the respective flux limits (see 
text). The adopted minimum luminosity of 2.55 X 10 44 /i^"q 2 ergs -1 
is indicated by the dashed line. Error bars are not shown. Note 
that even though no clusters above our adopted flux limit were 
found at redshifts 0.5 < z < 0.7, the lack of any such detections 
is significant in the analysis. 



our flux and luminosity thresholds for a given cosmology 
and mass-luminosity relation, we must be able to assign a 
probability to the luminosity observed from a cluster, given 
its true luminosity. Ultimately, the distribution is related to 
a Poisson distribution in the number of photons detected; 
however, the nontrivial conversion from photon count rate to 
unabsorbed flux and the variation in exposure times over the 
sky make the solution from first principles computationally 
difficult. Instead, we simplify the problem by assuming that 
the distribution P(L\L), where L and L are the observed and 
true luminosities, is normal with mean L. The variance is a 
function of the observed flux, which we estimate empirically 
for each sample by fitting the flux errors versus flux to a 
power-law model. We find power-law slopes of 0.52 for the 
BCS and REFLEX samples and 0.56 for MACS, consistent 
with approximately Poisson scaling. 



3.3 Other data 

The primary purpose of this paper is to present an anal- 
ysis based only on the X-ray luminosity function (XLF) 
data described above, along with the priors described in 
Section [4] However, we also show results obtained from 
other cosmological data for comparison purposes, along with 
results of a combined analysis, in Section [5] Specifically, 
the additional data consist of observations of galaxy clus- 
ter gas mass fractions (/ gas ), the cosmic microwave back- 
ground (CMB), and type la supernovae (SNIa) . The SNIa 
and / gaB results shown are identical to those in lAllen et al.l 
(2008). The / gas data are reported in that w ork and the 
SNIa results are derived from the compilation of lDavis et al.l 
l|2007h . whi ch includes data from the ESSENCE surve y 
(60 targets: IWood-Vasev et~ai1l2007i : iMiknaitis et al.ll2007r i. 
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the SNLS first year da ta (57 targets; lAstier et al-llioOr^ , 
45 nearby supernovae (| Jha et al.l 120071 ) and the 30 high- 
redshift supernovae di scovered by the H ubble Space Tele- 
scope and reported bv lRiess et al.l l|2007l ) for which a 'gold' 
rating was awarded (192 SNIa in total). Our analysis of the 
cosmic microwave background anisotropics uses three-year 
Wilkinson Microwave Anisotropy Probe (WMAP) data, in- 
cluding marginalization over a plausible range in the am- 
plitude of the S unyaev-Zel'dovich si gnal in galaxy clusters 
(0 < A S z < 2) <|Spergel et al.l 120071 ) . We use the October 
2006 version of the WMAP likelihood code0 

When analysing the CMB data, as with the XLF anal- 
ysis, we include the effects of density perturbations in the 
dark-energy fluid. For our standard analysis, we fix the 
dark-energy sound speed in units of the speed of light, c s , 
to 1.0. (Some discussion of the influence of perturbations 
can be found in Section 16.51 ) In addition, no prior on the 
scalar spectral index, n s , is used in the CMB analysis, since 
WMAP places very strong const raints on this parameter 
(0.951 ± 0.016; ISpergel et alT 2007); as noted in Section f6. II 
this degree of variation in n s produces negligible change in 
the XLF results. 



Table 1. Priors used in the analysis. All parameters not listed 
were assigned uniform priors on their physically allowed domains. 
N(fi, a) indicates the normal distribution with mean fi and vari- 
ance ex 2 , and U(a, b) indicates the uniform distribution with end- 
points a and b. 



Prior 


Parameter 


Density 


standard: 


w 


U(-5. 0,0.0) 




h 


N (0.72, 0.08) 




fl b h 2 


JV(0.0214, 0.002) 




n 3 


fixed at 0.95 




c 2 


fixed at 1.0 




A 


N (0.316, 0.0632) 




7 


U(—CI 35 35") 




'/i 


C/(-0.3,0.3) 




B 


7V(-0.125, 0.03) 






See Sectionl3.1.1l 


weak: 


h 


iV(0.72, 0.16) 




fl b h 2 


JV(0.0214, 0.004) 




n s 


iV(0.95, 0.05) 




c 2 


(7(0.0,1.0) 




.4 


JV(0.316, 0.1264) 




7 


C/(-0.7,0.7) 




</i 


f/(-0.6,0.6) 




B 


JV(-0.125, 0.06) 



4 ANALYSIS 

We parametrize the full model fitted to the X-ray luminosity 
function data as {h, Q, b h 2 , Q. c h 2 , as, n s , w, A, a, f3, 7, r?o, 
B}, where Q b and f2 c are the baryon and cold dark matter 
densities (S7 m = f2b+^c), A is the normalization of the Jenk- 
ins mass function (Equation [3} , and the last 6 parameters 
describe the mass-luminosity relation and its calibration. In 
addition to the assumption of sp atial flatness, we adop t the 
Gaussian priors h = 0. 72 ± 0.0 8 (Frccd man et al.ll200ll ) and 
VL b h 2 = 0.0214 ±0.002 l|Kirkman et alj|2003h from the Hub- 
ble Key Project and Big Bang nucleosynthesis studies, re- 
spectively. These latter priors are necessary because the like- 
lihood depends very weakly on h and £l b , which enter only 
through the transfer function. We additionally marginalize 
over a 20 per cent uncertainty in the normalization, A, of the 
J01 mass function to account for the residuals of the fitting 
formula to their simulations over the mass range of interest 
and the expected variation among cosmologies. We must also 
place a prior on the mass-luminosity evolution parameters 
7 and rfi, which are not constrained by our analysis of the 
RB02 data. For the standard set of allowances used in this 
paper, we adopt the uniform priors [ — y | < 0.35, which cor- 
responds to a limit of approximately 20 per cent evolution 
in the normalization of the mass-luminosity relation within 
redshift 0.7, and 1 771 1 < 0.3. Since the results are insensitive 
to the spectral index within a reaso nable range (Se ction l6.ll) . 
we fix n s = 0.95 in accordance with lSpergel et ail (|2007l ) for 
the standard analysis. The dark-energy equation of state was 
bounded by a uniform prior, —5 < w < 0. For the remaining 
model parameters, uniform priors on the physically allowed 
domains were used. These priors are summarized in Table [T] 
(labeled "standard" priors). The sensitivity of our results to 
the choice of priors is analyzed in Section 16.11 

Cosmological constraints were determined via Markov 



Chain Monte Carlo, employing the Metropolis algorithm; 
the calculations were performe d using a modi fi ed ve rsion 
of t he COSMOMC codjj o f iLewis fc Bridle] l|2002i . see 
also iRapetti et al.l [20051 . 12007I ). Multiple, randomly initial- 
ized Markov chains were produced for each set of results, 
and convergence to the posterior distribution was moni- 
tored using the Gelman-Rubin criterion, R, which mea- 
sures the ratio of betw een-chain to within-chain variances 
l|Gelman fc Rubinll 19921 ). as well as by visual inspection. Ac- 
ceptable convergence was defined by the requirement R— 1 < 
0.05. The log-likelihood of the data given a set of model pa- 
rameters is decomposed into the sum £ a p+£ v +£xlf, whose 
terms are described in the remainder of this Section. 



4.1 Mass— luminosity likelihood 

The parameters describing the normalization and slope of 
the mass-luminosity relation, a and f3, were constrained us- 
ing the x 2 statistic defined in Equation [5] Only the RB02 
data were used. As mentioned previously, we fix the value of 
5 to be 0.014, calculated in our reference cosmology using the 
method described in Section 13.1.11 This makes values of x 2 
from different steps of the Markov chain directly compara- 
ble, while retaining the de-biasing effect of the modification 
to \ 2 ■ (The reference cosmology with respect to which x 2 
is defined is unimportant, since the Markov chain is sensi- 
tive only to differences in the statistic.) The log- likelihood 
associated with x 2 is defined as £ a g — — x 2 /2. 

At each step of the chain, the estimated intrinsic dis- 
persion in the mass-luminosity relation, rfo 2 , is computed 
by Equation 1101 The contribution to the log-likelihood from 
the dispersion is the logarithm of the gamma distribution 
described in Section \3. 1.1 1 



http: / /lambda. gsfc.nasa.gov 



http: / /cosmologist.info/cosmomc/ 
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(k - 1) In (x) 



feln (0). 



(12) 



where k and 8 are the solutions to Equation [TT] for v = 84 
and F — 14, x — vrjo 2 /rjo, and where we have omitted a nor- 
malization term that does not depend on the model param- 
eters. This likelihood penalizes models for which the intrin- 
sic dispersion is far from the estimated dispersion measured 
from the RB02 data. 



4.2 Luminosity function likelihood 

The likelihood that N clusters with inferred luminosities in a 
range dL exist in a volume dV can in general be written as a 
Poisson probability plus a correction due to the clustering of 
halos with one another. Given that our sample covers a very 
wide survey area (~ 2/3 of the sky) and includes only the 
most luminous, and therefore rare, objects, the clustering 
term is negligible com pared to the pure Poisson term (e.g. 
iHu fc Kravtsovll2003h . If the plane of redshift and inferred 
luminosity is divided into non-overlapping cells, then the 
likelihood of our data is simply 



P(N 1 ,N2,...)=H 



(13) 



predicted in the jth cell, respectively. The log-likelihood is 
then 



E 



[Nj ln(Nj 



(14) 



where the model-independent term — ^\ ln(jV./!) has been 
dropped. 

If the cells are taken to be rectangular, with the jih cell 
given by zf 3 < z < z (2) and ^ L < tf\ then 



Above, P(L|A/5oo) is a log-normal distribution whose width 
is the intrinsic scatter in the mass-luminosity relation, rj(z), 
and P(L\L) is a normal distribution whose width as a func- 
tion of flux is modeled as a power law, as described in Sec- 
tion [331 When evaluating th e conversion between M 500 and 
M324m using the method of IHu fc Kravtsovl (|2003h ■ we as- 
sume a universal concent ration of 4 . at r adius r2oo, con- 
sistent with the results of iGao et all (|2007l ), who find that 
the concentrations of very massive clusters varies little with 
mass and redshift through z = 1. We have verified that in- 
cluding a log-normal scatter in concentration of ~ 0.15 (e.g. 
IComerford fc Nataraiar] 120071 ; ISchmidt fc Alien! \200H ) has 
no effect on our results. 

The sum over the second term in Equation 1141 reduces 
to the integrated number of predicted clusters in the detec- 
tion region of the surveys (i.e. within the redshift range and 
above the luminosity and flux thresholds), independent of 
binning. However, the first term in that equation is depen- 
dent on the choice of binning. To make optimal use of the 
data the bin size should be very small, in which case the in- 
tegrals in Equation[l5]can be approximated as the integrand 
multiplied by AzAL. By introducing a factor of d\jd\ = 1, 
where dL is the luminosity distance to the mean redshift of 
the bin, we can factor out the constant and cosmologically- 
invariant bin area Az(AL / d 2 L )\ the logarithm of this term in 
Equation Q3] is independent of model parameters and may 
thus be dropped, simplifying the log-likelihood to 



^ \ dz d L 



, , - dV dn „. 

+ / dzdL- r . (18) 

' dz dL 



where the summation is over detected clusters and the inte- 
gral extends over the detection region of the surveys. 



5 RESULTS 



Ni = 



J 2 ) 



dz 



dV(z) 
dz 



f dn(z, L) 

dij . 

dL 



(15) 



where V(z) is the comoving volume within redshift z. In the 
absence of intrinsic scatter in the mass-luminosity relation 
and measurement errors in the observed luminosities, the 
derivative of the comoving number density would be simply 



dn(z,L) _ dM 500 (L) dn(z, M 32 4m) rfM 324 „ 

/sky (Z, L)- 



dL 



dL 



dM s - 



dM 5 , 



(16) 



Here / S k y is the sky coverage fraction of the surveys as 
a function of redshift and inferred luminosity (i.e. flux), 
dn/dM-z24m is the Jenkins mass function (Equation |3} and 
M$oo(L) is the mass-luminosity relation (Equation [5j. The 
presence of scatter requires us to take into account that a 
cluster detected with inferred luminosity L could potentially 
have any true luminosity L and mass Msoo , with some asso- 
ciated probability. To calculate the predicted number den- 
sity correctly, we must therefore convolve with these proba- 
bility distributions: 



dh(z, L) 
dL 



= f sky (z,L) / dL P(L\L) 



(17) 



dM 500 P(L\M 5 



j dn(z, M 32 4m) rfM 32 4m 

dM 3 24m dMsoo 



5.1 ACDM constraints 

For a ACDM cosmology, we compare the joint Sl m -crg con- 
straints obtained from the BCS, REFLEX and MACS data 
sets individually, as well as their combination, in Figure [7] 
The marginalized constraints from the combination of the 
three cluster samples are f2 m = 0.28^'^ and erg = 0.78tg 
(Table [2). 

These constraints are in good agr eement with r ecent , 
independent results from the CMB (|Spergel et all l2007f) 
and cos mic shear, as measure d in the 100 Square Degree 
Survey dBeniamin et al.l |2007| ) (Figure [HJ and CFHTLS 



Wide field ( Fu et al. 20081 ). Our results are also in good 



overall agreement with previous findings based on the ob 
served X-r ay luminosity an d temperature functions of clus- 



ters ( e.g. lEke et all Il998l ; iDonahue fc Voitl Il999l; iHenryl 
|2000l: iBoreani et all I2OO1F ISeliakl 120021: lAllen et al 12003 ; 
iPierpaoli et all 120031 ; ISchuecker et al.l 120031 ; iHenrvl 120041 ). 

although the correction to the hydrostatic mass estimates 
employed in the present study leads to our result on ag 
being, typically, somewhat higher for a given value of Q m 
(see Section I6.3|l . Our result on Q m is in excellent agree- 
ment with curren t constraints based on cluster / gas data 
(|Allen et al.l 120081 and references therein) and the power 
spectrum of gala xies in the 2dF galaxy redshift survey 
l|Cole et al.l l2005h and Sloan Digital Sky Survey (SDSS) 
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Table 2. Best-fitting values and 68.3 per cent confidence intervals for the model parameters obtained from the luminosity function 
data. Our main results from this study occupy the first two lines; the remaining results are listed in the order that they appear in the 
text. a B=BCS, R=REFLEX, M=MACS. 6 Priors not specified below are standard (see Table [T). 1: all standard; 2: weak h and Cl b h 2 ; 
3: weak n s ; 4: weak A; 5: weak 7; 6: weak 771; 7: weak B; 8: b and Si, fixed at nominal values; 9: S and s;, fixed to 0; 10: fixed r]o = 0; 
11: survey luminosity measurement error= 0; 12: marginalized over c^; 13: no dark energy density perturbations. 
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Figure 7. Joint 68.3 and 95.4 per cent confidence constraints on 
Sl m and ct 8 for a ACDM model from MACS (red), BCS (blue), 
and REFLEX (green) individually, and their combination (pur- 
ple) using standard priors (Table[T]). Note that only the 95.4 per 
cent confidence regions are visible for the individual BCS and 
REFLEX data sets. 



Figure 8. Joint 68.3 and 95.4 per cent confidence constraints 
on n m and trg for a ACDM model using the combined X- 
ray luminosity function (XLF) data (purple) and our standard 
priors (Table [J. Also shown are in dependent constraints from 
the CMB (blue: ISpergel et al.|[2007h and cosmic shear (brown; 
iBeniamin et al.ll2007lV r 



(Eiso nstein et all 12003 ; iTegmark et alj 2006; Percival et al.l 
2007), as well as the com bination of CMB da ta with a variety 
of external constraints (|Spergel et alJl2007t ). Our result on 
as is marginally lower than that determined by weak lens- 
ing tomography in the Cosmic Evolution Survey (COSMOS; 



iMassev et al.l[2007l ) and by the observed number density of 
optic ally-selected groups and clusters in th e 2dF (|Eke et al.l 
l2006h and SDSS surveys ijRozo et alJl2007f ). 
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Figure 9. Joint 68.3 and 95.4 per cent confidence constraints on f2 m and erg (left) and f2 m and w (right) for a constant-ui model using 
the X-ray luminosi ty function data (pu rple) and sta ndard priors (Table [TJ. Also show n are independent constraints from the CMB (blue; 
ISpergel et al]|2007h . SNIa data (green; iDavis et al.|[2007t) and cluster / gas data fred lAllen et al.ll2008l) . and the combination of all four 
(gold). In the left panel, the dot-dashed lines indicate the XLF results using our weak prior on n B . 




Figure 10. Joint 68.3 and 95.4 per cent confidence constraints on f2 m and <rg (left) and erg and w (right) obtained from a combined 
/gas+CMB+SNIa analysis (blue) and the improved constraints obtained by combining these data with the XLF (gold). No priors on h, 
Q.^h? or n s are imposed in either analysis. In the left panel, the results from the XLF alone using standard priors (Table [TJ are shown 
(purple) in order to illustrate the degeneracy breaking. 



5.2 mCDM constraints 

5.2.1 Results using the XLF data alone 

The joint constraints on Q m and as from the luminosity 
function data using our standard priors (purple contours) 
are compared with those of WMAP (blue) in the left panel 
of Figure [5] Since the WMAP results are marginalized over 
n s , we also show the XLF results using our weak n s prior 
as dot-dashed lines, although the difference from the stan- 
dard results is small. The right panel displays constraints 



on fi m and w obtained independently from the XLF (pur- 
ple), WMAP (blue), SNIa data (green) and cluster / gas data 
(red) . Our new XLF results are consistent with each of these 
independent data sets, and with the cosmological-constant 
model (w = —1). The marginalized results from the X-ray 
luminosity function data are S7 m = ag = 0.851q'2q 

and w = -lA±° -$ (Table 01. 
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Table 3. Best-fitting values and marginalized 68.3 per 
cent confidence limits on cosmological parameters obtained 
from combined analyses of the / gas +WMAP+SNIa and 
XLF+/ g as+WMAP+SNIa data. No external priors on h, U h h 2 
and n s are used. 

/ g a S +WMAP XLF+/ gas 
+SNIa +WMAP+SNIa 

fi m 0.258 ± 0.022 0.269 ± 0.016 
<r 8 0.79 ± 0.06 0.82 ± 0.03 

w -0.99 ±0.07 -1.02 ±0.06 



5. 2. 2 Combined XLF+f sas + CMB+SNIa results 

We now consider the improveme nt over a com bined 
/ gas ±CMB±SNIa analysis (following lAllen et aLll200sT ) that 
can be achieved by additionally including the XLF data. The 
results of this comparison are displayed in Figure [10] and Ta- 
ble [3] The / gas ±CMB±SNIa combination already provides 
tight constraints on fi m , h, Q^h 2 and n s (hence no priors on 
these parameters are used in either combined analysis), but 
the degeneracy between w and erg (right panel of Figure [TIT)) 
limits the precision of the dark energy results. The addition 
of the XLF data breaks the degeneracy in the r2 m -er 8 plane 
(left panel), resulting in tighter constraints on f2 m , erg and 
w. 
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Figure 11. The number of objects N above detect flux S as a 
function of S (log N — log S distribution) for BCS (dotted line, 
above/right), REFLEX (solid line, above/right) and MACS (solid 
line, below/left). Shaded areas indicate the predictions and 68 per 
cent confidence regions of the best-fitting ACDM model. Due to 
the similar sky areas covered, the BCS and REFLEX predictions 
are essentially identical at S > 8. 



5.3 Goodness of fit 

In order to assess the goodness of fit for the XLF data, we 
compare the number of clusters in each sample with the 
number predicted by the best-fitting set of parameters in 
the Markov chains. Since the chains were generated at finite 
temperature, these do not necessarily correspond to the pos- 
terior modes, but they should be nearby. Using the XLF data 
only, the best-fitting ACDM parameters predict 73, 127 and 
37 clusters in BCS, REFLEX and MACS, respectively, while 
the best fitting mjCDM parameters predict 76, 130 and 35. 
Both are in good agreement with the true samples, which 
respectively contain 78, 130 and 34 clusters above flux and 
luminosity thresholds of Section T3. 21 This agreement for the 
ACDM case is displayed graphically as a function of the flux 
limit in Figure [TT1 (The numbers quoted above correspond 
to the low-flux ends of the lines in the figure.) The equivalent 
figure for wCDM looks similar. 

6 DISCUSSION 

6.1 Sensitivity to priors 

The standard priors used in our analysis are conservative; 
however, we have also performed analyses using the even 
more conservative priors labeled as "weak" in Table Q] These 
weak-prior results effectively demonstrate the degeneracies 
between parameters of interest (f2 m , erg and w) and each 
nuisance parameter, given the current level of statistical and 
systematic uncertainty. Two of these weak priors, on r\\ and 
B, did not result in significant changes to any of the cosmo- 
logical constraints. The results of the other tests are summa- 
rized below. One-dimensional, marginalized results for each 
of these analyses are also presented in Table [2] 
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Figure 12. Comparison of joint 68.3 and 95.4 per cent confidence 
constraints on n m and erg for a constant-ui model using standard 
priors (solid, black lines) and weak priors on h and Q^,h 2 (dashed, 
blue lines). Weakening these priors results in a slight expansion 
of the confidence region along the C m -crg degeneracy axis. The 
priors have no affect on the determination of w. 



6.1.1 Hubble constant and baryon density 

We consider h and Q.-oh 2 together, since both influence the 
results only through the shape of the power spectrum. The 
results obtained by doubling the width of the Gaussian pri- 
ors on these parameters (blue, dashed contours) is compared 
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Figure 13. Comparison of joint 68.3 and 95.4 per cent confidence 
constraints on £7 m and w for a constant-ui model using standard 
priors (solid, black lines) and a weak prior on the normalization of 
the mass function, A (dashed, blue lines). The constraints on Q m 
and w are weakened, although as is not greatly affected (Tabled- 



with our standard-prior constraint in Figure [T2j The effect 
of weakening these priors is primarily to expand the con- 
straints along the S7 m -crg degeneracy curve. The constraint 
on w is not affected. 



6.1.2 Scalar spectral index 

We next consider the effect of marginalizing over n s using 
a Gaussian prior of width 0.05, si gnificantly larger th an the 
constraint obtained from WMAP l|Spergel et al.l l2007). Like 
h and Qbh 2 , the effect of freeing this parameter is to weaken 
the constraints along the f2 m -crg degeneracy, although the 
effect is relatively small, even with this wide prior (left panel 
of Figure . 



6.1.3 Mass function normalization 

The weak prior on A doubles the uncertainty on the JOf 
normalization to 40 per cent. The result is a decrease in 
constraining power on S7 m and w, but not significantly on 
as - The expansion of the dark-energy constraints can be seen 
in Figure [TBI 



6.1.4 Evolution in the mass-luminosity normalization 

Since 7 describes only the redshift evolution of the mass- 
luminosity relation, doubling the size of its uniform prior has 
no effect on the determination of the redshift-zero quantities 
f2 m and as, as can be seen in Table [2] However, there is a 
degeneracy with w, resulting in a poorer constraint (Table[2] 
prior 5). 



6.2 AGN contamination 

The presence of active galactic nuclei (AGN) and other 
point-like X-ray emitters is a potential concern for the ac- 
curacy of the mass-luminosity relation. There is no pos- 
sibility of subtracting these point sources from the RASS 
data, since the number of photons in a cluster detection is 
typically too small. However, they are subtracted from the 
RB02 luminosities, to the extent that ROSAT can resolve 
them. In addition, there is the possibility of an increase in 
the density of AGN over the redshi ft range of our data (e.g. 
lHasinger. Mivaii. fc Schm idt 2005 and references therein). 
There may therefore be an inequivalence between luminosi- 
ties inferred from the RASS data and those used in the cali- 
bration of the mass-luminosity relation. In our work study- 
ing MACS clusters (in preparation), we take advantage of 
the Chandra X-ray Observatory's high spatial resolution to 
efficiently identify point sources and determine their con- 
tribution to the total cluster luminosity. We find that this 
contribution is at the few per cent level - much smaller than 
the typical survey flux uncertainty - and conclude that point 
source contamination is not an issue for the current study. 



6.3 Bias in mass measurements 

The bias in the measured masses of the mass-luminosity 
data is completely degenerate with the mass function nor- 
malization, A, which in turn affects inferences made on the 
other parameters. Consequently, it is important to include 
in the analysis the mean correction to masses motivated by 
departures from hydrostatic equilibrium, b, and the scatter 
in this bias among clusters, Sb, as well as the uncertainty 
in both quantities. Figure [T4l shows that the constraints ob- 
tained by fixing the mean bias and bias scatter at their nom- 
inal values (blue, dashed contours) are spuriously tighter 
than those of our standard analysis (black, solid contours) 
by ~ 50 per cent on £2 m and as and ~ 30 per cent on w. 
Also shown is the result of neglecting the bias correction en- 
tirely (red, solid contours), which results in lower values of 
as compared with the nominal correction. 

Comparing the results for priors 7 and 8 in Table[2]with 
our standard-prior results, it is clear that the uncertainty in 
Sj, has a greater effect than the uncertainty in b. This is 
encouraging, since the former is due primarily to the very 
small number of simulations used to calibrate the observa- 
tional mass bias. Straightforwardly increasing the number of 
simulated clusters should result in significant improvement 
to our constraints. 



6.4 Intrinsic and measurement luminosity scatter 

The mass distribution of clusters available to a flux-limited 
survey is influenced by the degree of scatter in flux for a 
given mass. As described in Section 14.21 this effect can be 
decomposed into convolutions due to intrinsic scatter in the 
mass-luminosity relation and measurement error in the sur- 
vey flux determinations. Failure to account for either of these 
sources of scatter when evaluating the number of detectable 
clusters predicted by a set of model parameters can signifi- 
cantly bias the result. The magnitude of the effect is demon- 
strated in Figure 1151 which shows the results for a raCDM 
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Figure 14. Comparison of joint 68.3 and 95.4 per cent confidence constraints for a constant-ui model using standard priors (black, solid 
lines), results with b and sj, fixed at their nominal values (blue, dashed lines), and results with no correction for hydrostatic mass bias 
(b = 0, s;, = 0; red, solid lines). The uncertainty in Sf, in particular has a significant impact on our results (see text). Neglecting the 
uncertainties on b and sj, produces spuriously tight constraints on Qm, erg and to. Neglecting completely the hydrostatic mass bias results 
in a significantly lower estimate of as- 
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Figure 15. Comparison of joint 68.3 and 95.4 per cent confidence constraints on f2 m and erg for a constant-ui model using standard 
priors (solid, black lines) and results obtained by ignoring either the intrinsic dispersion in the mass-luminosity relation (left) or the 
survey flux measurement errors (right) (dashed, blue lines). Failing to account for the intrinsic dispersion in the mass-luminosity relation 
produces a significant bias towards higher values of as and lower values of C m ; a similar, but smaller, bias is evident when the survey 
flux measurement error is not accounted for. 



cosmology obtained by ignoring either of these scatters in- 
dividually. In each case, the results are biased towards lower 
Q,m and higher erg, along the degeneracy between the two 
parameters. The constraints on cosmological parameters, in- 
cluding w, are also spuriously tight in both cases. (The joint 
Qm-w constraints are not shown, but closely resemble the 
tighter results in Figure [14]) This results from the fact that 
any scatter in luminosity compounds the systematic effect of 
uncertainty in b. In the case where intrinsic mass-luminosity 



dispersion is ignored, the contributions due to uncertainty 
in 770, 771 and Sb also vanish from the error budget. 



6.5 Dark-energy perturbations 

In the analysis thus far, we have accounted for dark-energy 
density perturbations assuming that the sound speed is 
equal to the speed of light (c^ = 1). Table [5] (prior 12) 
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shows that marginalizing over such a constant sound speed 
(0 < <? 3 < 1) has no effect on our results (but note that 
dark-energy models generally predict a time-varying sound 

speed). 

IWeller fc Lewis! (|2003h . iRapetti et all (|2005l ) and 
ISpergel et al.l l|2007l ) showed that including such perturba- 
tions greatly reduces the ability of WMAP data to discrim- 
inate among models with w < — 1 (Figure c/. Figures 
15 and 16 of Spergel et al.). In contrast, the effect of dark- 
energy perturbations on the XLF constraints is less dra- 
matic; Table [2] (prior 13) shows that an analysis with no 
perturbations results in somewhat larger constraints on w, 
opposite in character from the WMAP analysis. 



7 CONCLUSION 

We have presented new constraints on cosmological-constant 
(ACDM) and constant-u; (kiCDM) dark-energy models us- 
ing the observed X-ray luminosity function of the largest, 
most X-ray luminous galaxy clusters out to redshift 0.7, 
in combination with standard priors on h and ft^h 2 . At 
68.3 per cent confidence, we find fi m = 0.281q an d 
a 8 = 0.78±Q.i| for a ACDM model, and Q m = 0.24±°;^, 
as = 0.85±£^ and w = -1A±° ^ for a wCDM model. 
These results include marginalization over uncertainties in 
the theoretical mass function, non-self-similar evolution in 
the mass-luminosity relation, and the bias in mass estimates 
from X-ray observations. Our results constitute the first pre- 
cise determination of the dark-energy equation of state using 
measurements of the growth of cosmic structure observed 
in galaxy clusters, and provide additional support for the 
cosmological-constant model. 

These results build upon, and are largely in agreement 
with, a number of earlier galaxy cluster studies (see Sec- 
tion!]]). Our constraints on cosmological parameters are con- 
sistent with independent findings from studies of type la 
supernovae, anisotropics in the CMB, the X-ray gas mass 
fraction of galaxy clusters, galaxy redshift surveys and lead- 
ing cosmic shear surveys. The agreement between the results 
from these independent techniques is reassuring, and moti- 
vates a combined analysis of the data in order to investigate 
more complex models of dark energy. Combining our data 
with CMB, SNIa and / gas data, we find fi m = 0.269 ±0.016, 
a 8 = 0.82 ± 0.03 and w = -1.02 ± 0.06. 

The results for as presented here are somewhat higher 
than those from previous work based on the BCS and RE- 
FLEX data due to our correction to the masses used to 
constrain the mass-luminosity relation. The magnitude of 
this discrepancy underscores the need for an improved un- 
derstanding of the observational biases resulting from as- 
phericity, projection effects and hydrostatic disequilibrium. 
More advanced and comprehensive simulations, calibrated 
by gravitational lensing studies, show considerable promise 
in this area. Broad-band, high spectral resolution X-ray data 
from, for example, the New X-ray Telescope (NeXT) and 
Constellation-X, will allow precise measurements of gas ve- 
locities and non-thermal emission components in clusters, 
providing a more comprehensive understanding of the rel- 
evant gas physics. In the short term, simply increasing the 
number of simulated clusters used to calibrate the observa- 
tional bias would reduce the uncertainty in the scatter in 



bias among clusters. Since the uncertainty in this scatter is 
the dominant source of systematic uncertainty in as, this 
approach provides a relatively easy means of improving the 
precision of our results. 

With regard to dark energy, the systematics due to the 
theoretical mass function and the redshift evolution of the 
mass-observable relation are also important. The former can 
be addressed with a large suite of cosmological simulations, 
allowing a more precise parametrization of the cosmolog- 
ical dependence of the mass function. The latter necessi- 
tates a more rigorous study of galaxy cluster virial rela- 
tions and their evolution, carefully accounting for selection 
effects. Such improved analysis techniques will be required 
if future, high-redshift X-ray (e.g. Spectrum-RG/eROSITA) 
or Sunyaev-Zel'dovich surveys are to be used to their full 
potential. 
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